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Abstract 


Autogenous and allogenous control of the phytoplankton standing crop and 
biomass was proved in Mex Bay, west of Alexandria (Egypt). Among the clear signs of 
eutrophication are the heavy visible blooms, which became regular events, at 
intermittent periods during the warm seasons. The causative organisms progressed 
differently, the dinoflagellate, Gymnodinium catenatum, is well known toxic species. 

The Environmental Impact Assessment model (EIA) was applied to Mex Bay. It is 
a three-dimensional water quality and flow model, based on the three interacted 
dimensional: longitude, latitude and depth, calculations of time and space variations. 
The model deals with processes that affecting the water quality. It was initiated with 
measured field data and load input. The model was used to simulate phytoplankton 
biomass transport and ambient nutrient concentrations (nitrate, phosphate & ammonia). 
The model results and the field measurements are compared in order to fit the model. 
Generally, the calculated flow and concentrations corresponded well with observations. 
Deviations are attributed to several reasons. A notable influence of the main outfall on 
the concentrations, particularly in the near shore area, as well as, on the seawater quality 
was detected. The spatial distribution of the phytoplankton standing crop coincided 
clearly with the accompanied chlorophyll a concentrations calculated by the model. The 
results stressed the need to reduce load input into Mex Bay, in order to counteract 
eutrophication in Alexandria waters. Nitrogen reduction is more important for the whole 
area, while phosphorous reduction is more effective in the bay. 


Introduction 


Although coastal marine transiti - 
onal ecosystems located between land 
and open sea represent only small area 
of the world's ocean, they are of great 
ecological and economic importance. 

The land inflow of waste waters 
and subsequently the excessive input of 


ication. This phenomenon is of wide 
occurrence in various parts of the world 


(Vollenweider, 1981). Over 183x106 
m3 of partially treated or untreated 


sewage and waste waters are discharged 
annually into the coast of Alexandria, 


nutrients, mainly nitrogen and phosp - 
horus is underlying cause of eutroph - 
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but mostly into it's the western part. 
About 40% of the industry surrou - 
nds Alexandri city and as a conseque - 
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nce of the rapid increase of the human 
population and other activities, the task 
of controlling pollution have become a 
serious problem. The discharge waters 
rendered the neritic waters of 
Alexandria highly eutrophic, developing 
the occurrence of regular toxic and/or 
non toxic red tide blooms (Labib, 1992; 
Labib, 1994, a, b; Labib, 1996; Labib 
& Halim, 1995). Attempts have been 
made to evaluate the quantity of the 
different sources of pollutants, west of 
Alexandria (e.g. Said et al., 1991; Abdel 
Moneim et al., 1994; Nessim, 1994), 
and their ecological impacts on the 
standing crop (Dorgham et al., 1987; El- 
Sherif, 1989; Samaan et al., 1992). 


A case study: Mex Bay 


Mex Bay has a surface area of 
about 19.4 Km2, mean depth of about 
10 m and a volume of about 190.3x106 


m3(Labib, 1997). The bay receives 
directly daily agricultural discharge 


water, mean of 6.5x106 m3,(Emar „et.al 
1992) from neighbouring Lake 
Maryout, through Umum Drain. This 
water is mixed with municipal and 
industrial wastes. Five main industries, 
Chemical, Chloro-alkali, Tanneries, 
Cement and Petroleum, dispose their 
waste waters directly into the western 
part of the bay. Due to water 
circulation, it is also affected by 
additional amount of municipal 
wastewater from the main sewer of 
Alexandria located at its eastern part. 
Man-made eutrophication in Mex 


Bay caused drastic environmental 
changes. The bay has distinct structural 
dynamic properties that easily differe - 
ntiate it from the adjacent Mediter - 
ranean open oceanographic system. It is 
a site of heavy phytoplankton blooms. 

The aim of the present study is to 
set up the Environmental Impact 
Assessment (EIA) model to Mex Bay 
to: - provide an overview of the proce - 
sses governing the transport of the 
dissolved inorganic nutrients (ammonia, 
nitrate and phosphate); and the seasonal 
and spatial distribution of phytoplan - 
kton biomass (Chlorophyll a); - to 
determine the extent of eutrophication 
and to discuss the hydrodynamic 
characters of the bay. 

The results gained can give a 
better insight into practical processes 
have to be carried out in order to 
ameliorate the situation, if possible. 
The final goal of the work is not just to 
study of nutrient cycles and dynamics as 
well as biological processes, but to 
predict future trends of eutrophication 
in changed physical and chemical 
conditions. 

The bay is ideally suited for 
achieving the purposes of the study 
because its circulation patterns are a 
complicated function of winds, huge 
input of waste waters of different 
sources and pollutants and because of 
the ecological and economic problems 
associated with the progressive 
eutrophication. 

The study is based on surface data 
collected during 1992-1995. 


A Water Discharge (m3) 








" 2 م‎ 
e d o 
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Figure 1: water discharge from Umum Drain (Nessim,1994) 
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Nutrient concentrations and water 
discharge input 

The monthly amount of the 
discharged water (m°) into Mex Bay 
throughout Umum Drain was shown in 
Figure. 1. Data obtained from Ministry 
of Irrigation, Directorate of Mex 
Pumping Station, Alexandria. 

The daily nutrient injection from 
land-based sources led to continuous 
regeneration of nutrient. The high 
ammonia, nitrate and phosphate 
concentrations during May 1992 and 
February 1995 (ammonia reached 60.5 
u mol. lI! during the latter month), 
associated with peaks of discharge 
water input, and very limited 
phytoplankton growth. However, severe 
impoverishment occurred during the red 
tide bloom periods in July 1993 and 
October 1994 (ammonia and nitrate at 
4.2 and 3.5 u mol. 1! during the latter 
month, respectively). 


Phytoplankton standing crop, species 
composition and chlorophyll a 


The phytoplankton standing crop 
and species composition progressed 
differently during the different months. 
The density was relatively low in May 
1992, attaining 0.12.10° cell. I’. The 
community was mainly represented by 
microflagellate species, 56.54% to the 
total standing crop, followed by 
Scendesemus quadricauda (12.91%), 
Protoperidinium  steinii (4.9%) and 
Prorocentrum minimum (3.8%). 
However, chlorophyll a was relatively 
high (7.3 ug. 1^, influenced by the 
occurrence of fresh-water forms, 
mainly, Euglena acus and E. gracilis.. A 
water discoloration was observed on 12 
July 1993, the centric diatom, 
Skeletonema costatum, culminating a 
population size at 12.10° cell. I, 
96.2%.Chlorophyll a reached its 


maximum of 33.8ug. I". The accomp - 
anied minor species were Rhizosolenia 
fragilissima, Scrippsiella trochoidea and 
Prorocentrum triestinum. Again, 5-11 
October 1994 was a bloom period 
(maximum 44.2.10? cell.1), Gymno - 
dintum catenatum (53.13%), and 
Euglena granulata (26.56%), were the 
responsible species. Chlorophyll a 
content attained 27.3 ug. 11. These 
previous blooms covered the entire 
bay, but with different degrees. The 
western area was always dense. The 
phytoplankton standing crop was very 
poor in February (70.1.10° cell. I, 
despite the remarkable nutrients incr - 
ease, Phacus triquter contributed 51% to 
the total. Temperature was seen signi - 
ficantly influencing the density. 
Nutrients for their daily replenishment 
seen of negligible effect. 


The Environmental Impact Assess - 
ment (EIA) model 

Virtanen et al. (1994) describes 
the Environmental Impact Assessment 
(EIA) model in detail. 


The model calculations 
Water masses 


The water masses are treated as 
horizontal layers interfaced by levels at 
selected depths. In the horizontal plane 
the model area is subdivided into 
rectangles with arbitrary mesh intervals 
in both directions. Explicit finite- 
difference schemes are used for the 
numerical solution of flow velocities 
and of water level elevations. 
Mathematically, the physical interact - 
ions and mechanisms affecting water 
currents are described with Reynold's 
equations of motion, with the equation 
of hydrostatic pressure as a proper 
approximation of its vertical compon - 
ents and with the equation of 
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continuity(1): 


E =-V.p/ p, * V(KVju-2Q»u—u Vu F 


P-P +g. 1-2.0, T 1 p^ p,)Ewhere:u - 


rà 

di & dw 

ane 

flow velocity vector (u,v,w), t] -water 
level elevation, P -pressure, p, -average 
water density, Q -angular velocity of 
the earth rotation, p, -air pressure, K — 
dispersion matrix of momentum, g - 
gravity acceleration, F -external forcing 
vector = (0,0,-g). 

The water currents are calculated 
with the 3D multi-layer (Simons, 1980; 
Virtanen et al., 1986; Koponen et al., 
1992). 

The calculation of water currents 
is detached from that of transport and 
water quality (Virtanen and Koponen, 
1985; Halfon et al., 1990; Koponen et 
al., 1992). 


The flow fields UL induced by 
winds W and water discharges Q are 
calculated in advance from several 
levels of water surface 7. When the 


water level is changed, the flow fields 
for transport calculation changed 
accordingly. The surface area covered 
with water is updated. For any water 
level 77 between two levels of flow field 
calculatiz7, and7,,, the flow field U- 
ie. the flow velocities u = (u,v,w) at 
every location (m,n,k) of the 3D grid - 
will be : 


U= 1— Tha U+ 7), — 7] y+” (2) 


Ni Mh Ni — ha 
Calculation of water quality 

Release, transport, mixing, 
resettling and the changes of water 
quality and bottom properties are 
coupled with each other in a most 
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complex net of interactions. In the 
model these are gathered within the 
transport equation, with internal loading 
as one of its boundary conditions (Eq. 
3): 








d a gd a A æ Aa) Aa 
a la Yg DDE 


+A gl)! d-R Tc...] 

GL,=E,+ E+E, 

Where 

c- Concentration (or strength) of any 
quantity computed, yy sọ -flow 
velocity components obtained from Eq. 
23. D, D. -dispersion coefficient 
components of concentrations, qL - 
loading discharges (from point source P, 
rivers R, drainage area D, bottom B and 
through water surface A), n - discharge 
direction co-ordinate, R - changes from 
internal (biogeochemical) reactions 
within water. 

In the solution, the simultaneous 
complex interactions are governed by 
considering each process after the other 
with its immediate consequences and 
going through all of the processes step 
by step with short intervals until the 
desired time period as the sum of the 
time-steps for each process is equally 
expired (Virtanen et al., 1994). 


Model data and initial conditions 

The in-flowing Mediterranean 
waters at the boundaries were calculated 
with an average current velocity of 8 cm 





secl. The wind data was obtained from 
the Egyptian Navy Meteorological 
station for the different months (Fig. 2). 


The model application in Mex Bay 

The model was calibrated for 4 
different seasons (May 1992, July 1993, 
October 1994 and February 1995), on a 
point, where measurements were 
available. Modelled results and field 
data were compared. 
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Results and Discussion 

For numerical simulation Mex 
Bay is subdivided with equal mesh 
resolution of 300 m into 40 by 17 
horizontal grid squares and in the 
vertical direction into 10 one meter 
thick layers. The amount of the 
discharge water from the Umum Drain 
is given as monthly constant input data 
for the model. The open boundary to the 
Mediterranean has a monthly constant 
concentration. The model bathymetry 
was obtained from the Admiralty Chart, 
by linear interpolation of the depth data 
into the model grid. 

The simulation results of the 
nutrients and chlorophyll a are 
presented as pictures of horizontal 
concentration distribution on surface. 
Wind frequency distributions are shown 
in Fig. 2. 

During May 1992, the 
predominant wind direction was NW 
and relatively less NE. In July 1993, the 
winds were less variable, mostly NW. A 
noticeable change occurred in October 
1994, when the wind direction was 
dominated by NE, relatively less 
powerful in NW. Again in winter 
(February 1995), the strong winds were 
mostly of NW direction, less affective 
towards south and east. 

The flow velocities, induced by 
several combination of wind, jet 
discharge water from Umum Drain and 
discharge water at the western part of 
the bay, as well as the partial effect of 
the main sewer of Alexandria at its 
eastern vicinity, are calculated with a 
time-step of 30s for periods of almost 2 
weeks. 

To show the differences in the 
resultant flow patterns under the wind 
action, a case study is presented. The 
water mass is subjected to the free 
surface shear induced by a uniform and 
steady wind of 5-m sec", blowing from 
the main eight directions. The steady 


state circulation patterns at the sea 
surface are sketched in Fig. 3. 

The simulation results of the 
model run for the different months are 
shown in Figures 4-8. The main 
indicators of water quality in the present 
study were the nutrients (ammonia, 
nitrate and phosphate). The phytopl - 
ankton standing crop was expressed as 
biomass. 

The results are shown as time- 
series of their concentrations, as 
concentration maps, and as numerical 
values, ranging from smallest to 
maximum on the computer screen, 
expressing the development of 
concentrations with colours along a 
selected horizontal lines. 

The spatial distribution of the 
nutrients and chlorophyll a obtained 
from the model results indicated that: 

e under the north-westerly wind in May, 
the maximum concentrations were 
restricted to the area around the outlet. 
However, very high nitrate and 
ammonia concentrations were also 
seen to spread towards the west. The 
time series at the selected points 
reflected such pattern. 

the predominance of NE wind 
direction in July affected clearly the 
distribution patterns, maximum 
concentrations pushed to a narrow 
coastal zone, and a gradual decrease 
far from the shore and towards west. 
The time series data showed 
variability in ammonia and nitrate 
concentrations, while phosphate 
exhibited a very limited range of 
variation. 

in October, following the prevailing 
wind directions, three distribution 
patterns were observed. Variations of 
phosphate with the time series were 
still limited 

the jet effect of land run-off on 
ammonia distribution was seen in 14 
February under the SN wind direction. 
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High concentrations covered most of 
the bay. The distribution patterns 
changed on 17 February, with the 
reverse wind direction, maximum 
concentrations were located near the 
coast, extending towards the inner 
part of the bay. 

e the characteristics of nutrient distrib - 
utions on the surface were in signif - 
icant connection with wind speed and 
direction. 

* the spatial distribution of chlorophyll 
a was clearly influenced by the trans - 
port processes and the jet effect of the 
discharged water. High concentrations 
were detected in front of the main 
outlet, as well as towards the western 
part of the bay, decreasing towards the 
open sea. 

As shown in Fig. 9, the correspo - 
nding spatial distribution of the standing 
crop was somewhat similar to the 
modelling results of chlorophyll a. 
However, noticeable deviations occur, 





May 1992 


N 





July 1993 


mainly attributed to the effect of 
changing species composition, excha - 
nge with the open sea flora and that 
coming from Lake Maryout. A compa - 
rison of the model results and field data 
is shown in Fig. 10. In general term, the 
behaviour of the model results is rather 
similar to the field observations, and the 
calculated values correspond well with 
the real data, but their values were 
lower, and there are some differences 
between them. These might be for some 
reasons; *the model results represent 
the mean values of the grid squares of 
the whole bay. 

*local irregulations in sampling can’t be 
avoided. 

*other land-based sources of nutrients 
in the bay are not included. 

*water exchange at the model open 
boundaries, significantly influence the 
concentration values in the modelled 
area. 





February 95 


Figure 2. Wind frequency (%) 
disrtibution over Mex Bay 








Figure 3. Modelled current field under 
a uniform and steady 5m/sec wind from NW 
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*vertical mixing in the model is more 
intense compared to nature. 
Subsequently, the concentrations are 
distributed more rapidly in the vertical 
direction, causing reduction of concent - 
rations in the surface layer. 

*The ecological factors affecting the 
concentrations are not taken into 
account. 

*the high supply from Umum Drain 
determines a prevalent SN gradient and 





a weaker NW gradient, persisting most 
of the time, but modified at times by 
prevailing meteo-climatic conditions. 
This powerful outflow creates a strong 
offshore surface current, and as a result 
there is a subsurface compensation 
current of higher salinity towards the 
coast (Labib, in preparation). 

the species composition which creates 
sharp deviations between the counts and 
biomass. 
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Figure 4. Model results. Horizontal distribution and time-series 
during May. 
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Figure 5. Model results. Horizontal distribution and time-series 
during July 
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Figure 6. Model results. Horizontal distribution and time-series 





during October. 





Transport Processes In A Heavy Eutrophic....... 







EJ EI a2 as De 


© 10 2030 40 50 60 70 80 90 100 








uM 


OUTPUT DEPTH: 
SCALE 





Figure 7. Model results. Horizontal distribution and time-series 
during February 1995. 














Figure 8. Model results. Horizontal distribution and time-series 
of chlorophyll a. 
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Figure 9. Spatial distribution of the standing crop 
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Figure 10. Modelling and field data measurements. 
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Conclusions 


The hydrodynamic and water 
quality model (the Environmental 
Impact Assessment model, EIA), with 
available field data, was applied to the 
highly eutrophic marine basin (Mex 
Bay). The bay, subjecting to a huge 
volume of discharge water, has distinct 
physical, chemical and biological 
aspects that easily differentiate it from 
the adjacent Mediterranean waters. The 
bay is a site of heavy algal blooms. The 
model was run for 4 different months 
during the period from 1992-1995, in 
order to follow the patterns of the 
spatial distribution of the dissolved 
inorganic nutrients, ammonia, nitrate 
and phosphate, and chlorophyll a. 
Generally, the calculated flow and 
concentrations corresponded well to 
observations and the model reacts 
logically to environmental changes. The 
model describes the seasonal variations 
of nutrients and chlorophyll a rather 
well. The largest deviations in the algal 
biomass results are observed in summer 
and autumn, affected mainly by the 
variability in species composition. A 
notable influence of the main outfall on 
the concentrations, reflecting a jet 
effect, was distinct. The concentrations 
decreased away from shore and towards 
the open sea. The western part of the 
bay, affected by transport processes, 
was denser and richer. Deviations are 
attributed to several reasons. Finer 
vertical resolution or stratification 
would be taken into account to reduce 
the vertical mixing. In future, the 
model must be tested and verified with 
extensive nutrient and phytoplankton 
data. 

It is recommended that the load 
reduction of nutrients, particularly 
nitrogen and phosphorous, and the 
water discharge input are very important 
to overcome serious health, economic 


and environmental problems, not only 
in Mex Bay, but also along the coast of 
Alexandria. The increased frequency of 
visible algal blooms, including toxic 
species, associated with heavier Man- 
made eutrophication could lead to 
undesirable situations. 
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تقييم عملية زيادة خصوبة مياه خليج المكس بالإسكندرية 
باستخدام نمودج تجريبي للتقييم البيئى 


*د/ سامية كامل ميخائيل- * a)‏ / وجدى لبيب جرجس- **د/إيمان 


إبراهيم صيام 
*معهد القومى لعلوم البحار 
** الاكاديمية العربية للعلوم والتكنولوجيا والنقل البحرى قسم حماية البيئة البحرية 


تناقش هذه الورقة البحثية موضوع زيادة خصوبة مياه خليج المكس قرب 
الإسكندرية وظهور ذلك على فترات متقاربة في فصلى الربيع والصيف وظهور 
السوطيات السامة التي تسبب كثيرا من التلوث. 

تم مناقشة هذا الموضوع من خلال نموذج حسابي لتقييم التأثيرات البيئية معتمدا 
على الأبعاد الأساسية الثلاثة الطول والعرض والعمق وكذلك على العمليات التى تؤثر 
على حالة مياه البحر. 1 

ويبدأ بقياس العوامل البيئية ثم تأثير العوامل الخارجيةء ويستخدم كذلك لتحفيز 
الكتلة الحيوية للهائمات النباتية وحركتها وكذلك تركيز نسبة الأملاح المغذية من 
أمونيا ونترات وفوسفات ثم مقارنة النتائج الحقلية والنتائج التي حصلنا عليها من 
النموذج التجريبي لنصل إلى النموذج الأمثل. 

وجاءت النتائج لتثبت أن القياسات الحقلية كانت ملائمة تماما مع الدراسات 
والملاحظات أما الانحرافات في بعض النتائج فكانت رد فعل للتغيرات المناخية نتيجة 
لاختلاف الفصول . 
وأظهرت النتائج أيضا أن توزيع الهائمات النباتية كان ملازما لتركيز الكلوروفيل 
(Nog oa‏ المحسوب بواسطة النموذج القياسي. 

وأوضحت النتائج ضرورة ة السيطرة ة على الملوثات الخارجية Gx uad‏ فين 
خليج المكس لتقليل حالة زيادة الخصوبة وذلك عن طريق تقليل زيادة تركيز 
النيتروجين والفسفور في مياه الخليج. 
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